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ABSTRACT 

It is shown that differential Doppler shift of different patches of the blastwave 
front in the fireball model at varying angles to the line of sight could provide 
an intrinsic smoothing mechanism for the spectra of gamma-ray bursts (GRBs) 
and the associated afterglows at lower energy bands. For the model parameters 
of interest, it is illustrated that a monochromatic spectrum at v in blastwave 
comoving frame is smoothed and observed to have a half-width-half-maximum 
(HWHM) of ~ (0.6 — 1SS)v. Some other imphcations of this smoothing process 
are discussed. In particular, if the circumburster medium is uniform and electron 
and magnetic energies are fixed fractions of the total post-shock energy with 
time, the observed GRB or afterglow spectra cannot be steeper than with 
a = 0.75 — 1.25, regardless of the intrinsic spectra in the comoving shock frame; 
i.e., a very steep electron distribution function (with p > 3) could produced an 
observed spectrum with a ~ 1. In addition, a generic fast- rise-slow-decay type 
of GRB temporal profile is expected. 

Subject headings: gamma rays: bursts - hydrodynamics - relativity - shock 

waves 
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1. Introduction 

It is remarkable that the simple fireball model for cosmological gamma-ray bursts can 
explain the many major features of the gamma-ray bursts and their afterglows (Rees & 
Mezsaros 1992, 1994; Paczynski & Rhodes 1993; Wijers, Rees, & Mezsaros 1997; Vietri 
1997a,b; Waxman 1997a,b; Reichart 1997; Katz & Piran 1997; Sari 1997). It is noted, 
however, that the fireball model may produce multiple spectral components due to the 
existence of both the forward and reverse shocks at least at some time period (e.g., 
Mezsaros, Rees, & Papathanassiou 1994). Moreover, reprocessing of the primary spectrum 
by, for example, inverse Compton scattering (e.g., Pilla & Loeb 1998), introduces additional 
features to the spectra, even if it is initially featureless. Finally, resonant line features or 
recombination lines due to heavy elements, for example, in dense blobs as discussed in 
Mezsaros & Rees (1998), may add additional features to the continuum. In this Letter it is 
pointed out that spectral smoothing due to differentially varying Doppler shift for different 
patches of the fireball at varying angles to the line of sight provides an intrinsic, unavoidable 
smoothing mechanism. While this effect is fairly well known, a more quantitative analysis 
focusing on the case of GRBs is useful. This Letter presents such a semi-quantitative 
analysis. 



2. Differential Doppler Shift Across a Fireball Front 

For the present illustration I assume that the GRB fireball is spherical and 
homogeneous, for which a single parameter, 9, the angle to the burster-observer vector, 
is sufficient to characterize the direction of a traveling patch of shock heated, radiation 
emitting material. It is convenient to use the time measured in the rest-frame of the burster, 
t, as the independent variable to express other quantities. 
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First, one has to find the relation between time t for a fireball patch at 9 and the time 
measured by the observer on the Earth (called "O" hereafter), tots, i-e., the arrival time. 
Since the apparent perpendicular traveling speed of a patch with 9 seen by is 

_ f](t) sing 
^^^^^ - 1- Pit) cos 9 

in units of the speed of light, where f3{t) is the spherical expansion speed of the fireball in 
the rest-frame of the burster. By definition, P±{t) is 

/3^{t)^-^sm9. (2) 



Combining equations (1, 2) gives 

dr (5{t) 



dtobs l — /3{t)cos9' 
Since one also has the following relation 



(3) 



one finds the equation relating tobs to t: 

dt 



flv 



= 1 -/3(i)cosg. (5) 



Next, to have a tractable treatment, it is assumed that the Lorentz factor, 
T{t) = IZ-^I — (3{ty, has the following simplified evolution: it is constant (equal to Fj) at 
t < tdec and decays at t > tdec as 

r = r,(-^)-", (6) 

''dec 

where Fj is the initial Lorentz factor and tdec (measured in the rest-frame of the burster) 
characterizes the transition time after which deceleration of the fireball expansion becomes 
significant (hence the fireball kinetic energy can be converted into radiation) and can be 
expressed approximately as t^ec — ( ^^pg^f^ n )^^^ (Blandford & McKee 1976), where E is 
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the initial fireball energy, n is the density of the circumburster medium (which is assumed 
to be uniform, for simplicity) and other notations are conventional. Note that a = 3, if the 
fireball cools radiatively efficiently, and a = 3/2, if the fireball cools only adiabatically. 

To make the point in a simple way it is assumed that the total luminosity per unit 
frequency (in the comoving frame) of the shock front at a fixed time t is a delta function in 
frequency (i.e., a monochromatic spectrum): 

U{u',t)^C{t)S[u'-uo{t)] (7) 

in the frame comoving with the blastwave. The characteristic frequency i^o(^) at time t in 
the comoving frame is parameterized as 

Mt) = A{^r, (8) 

where A is a constant. Note that, in the case of synchrotron radiation and assuming that 
both the electron thermal energy density and the magnetic energy density are fixed fractions 
of the post-shock nucleon thermal energy density, one has ip = 3. C{t) is expressed as 

C{t) = B{-L)i, (9) 

''dec 

where B is another constant and the ^ parameterizes the temporal profile of the amplitude 
of the radiation. 

At a given time receives radiation from different parts across the fireball surface, 
emitted at varying times in the burster frame; i.e., radiation from regions of varying 9 at 
varying t [see equation (5) for the relation between t and tobs] is seen by O at the same time 
(Sari 1998). The received frequency of the radiation at tobs from region with 9 emitted at t 
is 

iy{tobs)^Mt)D{9,t), (10) 
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where D{9, t) is the Doppler factor for regions with 9 at time t: 



The flux density observed by O at frequency u at time tobs is 

S{iy, tobs) = ^ /_'^ t)D\e, t)dii (12) 

where v/D) is the frequency in the blastwave frame, d is the distance of the GRB from 
O and ji = cos^. Combining equations (10,11) gives 

Inserting equations (7,13) into equation (12) and integrating over u give 

^^"'"'^''^ - 8^^d?W) ^ ' ^ — m — ■ ^ ^ 

Note that S{i', tobs) is in a parametric form; given tgbs, one can determine the burster frame 
time t for a given patch at 9 using equation (5). Then, one determines 5" (z/, tobs) using 
equation (14) combined with equations (6,8,9,11), given t and 9. Meantime, i^(to6s) is 
related to t and 9 by equation (10). Thus, one can find S{i',tobs) as a function of i'{tobs) at 
tobs- Let us consider a few simple but relevant cases to illustrate the effect. 

Case 1: q; = 3/2,'0 = O and ^ = 0. This case may have some bearing on such radiation 
features as atomic line features whose intrinsic frequencies are independent of t (i.e., = 0). 

Case 2: a — 3, t/j — 3 and ^ = 1. This case may be related to the epoch where radiative 
cooling is efficient and intensity at the frequency in question is still rising, which could 
be relevant for radiation at frequencies lower than the peak of the synchrotron radiation 
spectrum (due to a truncated power-law electron distribution) at early times of a GRB 
event. 

Case 3: o; = 3, = 3 and ^ = —1. This may be related to the epoch where radiative cooling 
is efficient and intensity at the frequency in question has started to decrease, which could 
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be relevant for radiation at frequencies higher than the peak frequency of the spectrum at 
early times of a GRB event. 

Case A: a — 3/2, ip — 3 and ^ = 1, This case is similar to Case (2) with the primary 
difference that radiative cooling is unimportant here. This may be relevant for GRB 
afterglows such as the radio afterglows when electron cooling time is likely to be significantly 
longer than the dynamic time of the expanding fireball at frequencies lower than the peak 
frequency of the spectrum. 

Case 5: a — 3/2, t/j — 3 and ^ = — 1, This case is similar to Case 4 but for frequencies 
higher than the peak frequency of the spectrum. 

While it is convenient to express various quantities using t as the independent time 
variable, one needs to express the final observables using tgbs, which is related to t (for 
9 = 0; see equation 5) as 



for the simphfied solution of r{t) given by equation (6). Figure (la) shows the flux density 
as a function of frequency at tabs — tdec/^'i for the five cases. The frequency is normalized 
such that unity corresponds to the radiation from regions with = and the flux density is 
normalized to be unity at unity frequency. The sharp turns to the left for cases (ii,iii,iv,v) 
correspond to the sharp turn of the evolution of the Lorentz factor at i = tdec- Figure 
(lb) shows the flux density as a function of frequency at tobs = 5000idec/r^ for the flve 
cases. Also shown in both panels are two straight lines in the upper right corner indicating 
the spectral slope of —0.75 and 1.25, respectively, which bracket the range of the spectral 
slope for various cases shown. Note that sharp turns as seen in (la) are not visible simply 
because they appear at much lower intensity level than the displayed range in the figure. 



and n, tobs is of order a second and a day, respectively, after the fireball explosion for the 




(15) 




1/3 



days. Therefore, for typical values of Fj 
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two cases shown in (la) and (lb). These two cases may respectively be relevant for bursts 
in gamma-ray and afterglows at lower energy bands. It should be noted that, although the 
external shock model is used to illustrate the smoothing magnitude, the results should be 
applicable to the internal shock model as well. 

Recall that in all cases a delta function spectrum at v is assumed in the comoving 
frame at a given time. One sees that this delta function spectrum is smoothed out to appear 
as a broad spectrum with a half-width-half-maximum (HWHM) of ~ (0.6 — 1SS)v for all 
the cases considered, except for Case 1 where HWHM is ~ O.lz/. An immediate implication 
from Figure 1 (see the dot-long-dashed curves on the upper right corner in the two panels) 
is that the observed spectra of GRBs or their afterglows should roughly have ^ if the 
electron distribution function power index p is equal to greater than 3. In other words, the 
observed spectra cannot be steeper than regardless the value of p. This slope se 

It is of interest to understand where the different radiation that correspond to different 
frequencies in Figure 1 comes from. Figure 2 shows the frequency seen by O at a given time 
as a function of the emission shell radius from the burster, r, expressed in units of ctdec- It 
is seen that, for realistic cases (ii,iii,iv,v) that correspond to the continuum radiation of the 
GRBs and afterglows, higher frequency radiation comes from earlier time t (in the burster 
frame) with larger angle 9 up to tdec after which there is a downturn to lower frequency at 
still earlier times due to assumed constancy of F thus constancy of uo{t)- For case (i) lower 
frequency radiation comes from earlier time t with larger 6. In both panels (a,b) also shown 
along the solid curves using solid dots are the corresponding 9 values in degrees. Dotted 
and dashed curves are also punctuated by open circles and open squares, corresponding to 
the same 9 values. 
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3. Discussion 

It is shown that differentially varying Doppler boost of different patches of the fireball 
front provides an intrinsic, unavoidable smoothing mechanism for the spectra of gamma-ray 
bursts and their afterglows. The detailed smoothing patterns are complicated, depending 
upon various factors such as the evolution of the Lorentz factor and the evolution of the 
intrinsic (i.e., comoving frame) radiation spectrum. Nonetheless, for plausible ranges of 
model parameters of interest, a comoving frame delta function spectrum at v is smoothed 
to have a HWHM of ~ (0.6 — 2.0)i/, assuming that the time evolution of the characteristic 
frequency v{i) is proportional to some positive power of the shock front Lorentz factor, 
(where > 0; see equation 8). This type of smoothing may be applicable to continuous 
spectra such as from synchrotron mechanism. 

In the case = (appropriate for atomic line features which are independent of the 
blastwave dynamics), the spectral smoothing is smaller, with a HWHM of ~ ^.\v. Thus, a 
sharp linelike emission feature (in the comoving frame) would be smoothed out to have an 
equivalent width of about ^.\v. Furthermore, the spectral profile of such a feature will be 
asymmetrical with a sharp cutoff at the high end (see the two solid curves in Figure 1). 

Two interesting and natural consequences arise due to the differential Doppler 
smoothing. First, the observed spectra of GRBs or their afterglows cannot be steeper than 
v~'^ with a — 0.75 — 1.25 (Figure 1), even though the intrinsic spectra in the comoving 
shock frame may be much steeper, if the circumburster medium is uniform and electron 
and magnetic energies are fixed fractions of the total post-shock energy with time. Second, 
a generic fast-rise-slow-decay type temporal profile of the GRB bursts is expected, not 
necessarily reflecting the intrinsic temporal profiles of the bursts in the comoving frame. 
This can be easily seen by considering the case where the intrinsic (blastwave frame) 
spectrum is a bivariate delta function in both time and frequency. In this case the fast rise 
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occurs when the radiation from the region around ^ = enters the observer's finite band. 
Subsequently, radiation from regions with gradually increasing 6 is received at decreasing 
amplitudes in the same observer's finite band with a decaying time scale of ~ tdecf^^^ (in 
observer's frame). 

The work is supported in part by grants AST9318185 and ASC9740300. I thank 
Bohdan Paczyhski for discussion. 
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Fig. 1. — Panels (a,b) shows the flux density as a function of frequency for the five cases 
(see text), at tots = tdecl^f and tofes = SOOOt^ec/r?! respectively. A fiducial value of Fj = 300 
is used. In both panels there are two straight lines in the upper right corner showing the 
spectral slope of —0.75 and 1.25, respectively, which approximately bracket the range of the 
spectral slope for various cases shown. 

Fig. 2. — shows the frequency seen by (9 at a given time as a function of the emission shell 
radius from the burster r, expressed in units of ctdec for the five cases, with panel (a) at 
iohs = idecl^l and panel (b) at tobs = SOOOtdec/r? • A fiducial value of Fj = 300 is used. In 
both panels also shown along the solid curves using solid dots are the corresponding 9 values 
in degrees. Dotted and dashed curves are also punctuated by open circles and open squares, 
corresponding to the same 6 values. 
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